Adiabatic quantization of Andreev levels 
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We identify the time T between Andreev reflections as a classical adiabatic invariant in a ballistic 
chaotic cavity (Lyapunov exponent A), coupled to a superconductor by an iV-mode constriction. 
Quantization of the adiabatically invariant torus in phase space gives a discrete set of periods T n , 
which in turn generate a ladder of excited states e nm = (m + l/2)nh/T n . The largest quantized 
period is the Ehrenfest time To = A -1 In N. Projection of the invariant torus onto the coordinate 
plane shows that the wave functions inside the cavity are squeezed to a transverse dimension W/vN, 
much below the width W of the constriction. 
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The notion that quantized energy levels may be as- 
sociated with classical adiabatic invariants goes back to 
Ehrenfest and the birth of quantum mechanics 0. It 
was successful in providing a semiclassical quantization 
scheme for special integrable dynamical systems, but 
failed to describe the generic nonintegrable case. Adi- 
abatic invariants play an interesting but minor role in 
the quantization of chaotic systems 0, 0] ■ 

Since the existence of an adiabatic invariant is the ex- 
ception rather than the rule, the emergence of a new one 
quite often teaches us something useful about the sys- 
tem. An example from condensed matter physics is the 
quantum Hall effect, in which the semiclassical theory is 
based on two adiabatic invariants: The flux through a 
cyclotron orbit and the flux enclosed by the orbit cen- 
ter as it slowly drifts along an equipotential |4j. The 
strong magnetic held suppresses chaotic dynamics in a 
smooth potential landscape, rendering the motion quasi- 
integrable. 

Some time ago it was realized that Andreev reflection 
has a similar effect on the chaotic motion in an electron 
billiard coupled to a superconductor [j| . An electron tra- 
jectory is retraced by the hole that is produced upon ab- 
sorption of a Cooper pair by the superconductor. At the 
Fermi energy Ep the dynamics of the hole is precisely the 
time reverse of the electron dynamics, so that the motion 
is strictly periodic. The period from electron to hole and 
back to electron is twice the time T between Andreev 
reflections. For finite excitation energy e the electron (at 
energy Ep + e) and the hole (at energy Ep — e) follow 
slightly different trajectories, so the orbit does not quite 
close and drifts around in phase space. This drift has 
been studied in a variety of contexts 0, IS IE 13 1 but 
not in connection with adiabatic invariants and the asso- 
ciated quantization conditions. It is the purpose of this 
paper to make that connection and point out a striking 
physical consequence: The wave functions of Andreev 
levels fill the cavity in a highly nonuniform "squeezed" 
way, which has no counterpart in normal state chaotic or 
regular billiards. In particular the squeezing is distinct 
from periodic orbit scarring ^(j an d entirely different 
from the random superposition of plane waves expected 



for a fully chaotic billiard |Tl| . 

Adiabatic quantization breaks down near the excita- 
tion gap, and we will argue that random-matrix the- 
ory can be used to quantize the lowest-lying exci- 
tations above the gap. This will lead us to a formula 
for the gap that crosses over from the Thouless energy 
to the inverse Ehrenfest time as the number of modes in 
the point contact is increased. 

To illustrate the problem we represent in Figs. ^ an d 
|5] the quasiperiodic motion in a particular Andreev bil- 
liard. (It is similar to a Sinai billiard, but has a smooth 
potential V in the interior to favor adiabaticity.) Fig. Q] 
shows a trajectory in real space while Fig. [21 is a section 
of phase space at the interface with the superconductor 
(y = 0). The tangential component p x of the electron 
momentum is plotted as a function of the coordinate x 
along the interface. Each point in this Poincare map cor- 
responds to one collision of an electron with the interface. 
(The collisions of holes are not plotted.) The electron is 
retroreflected as a hole with the same p x . At e = the 
component p y is also the same, and so the hole retraces 
the path of the electron (the hole velocity being opposite 
to its momentum). At non-zero e the retroreflection oc- 
curs with a slight change in p y , because of the difference 
2e in the kinetic energy of electrons and holes. The re- 
sulting slow drift of the periodic trajectory traces out a 
contour in the surface of section. The adiabatic invari- 
ant is the function of x,p x that is constant on the con- 
tour. We have found numerically that the drift follows 
isochronous contours Ct of constant time T{x,p x ) be- 
tween Andreev reflections [l^- Let us now demonstrate 
analytically that T is an adiabatic invariant. 

We consider the Poincare map Ct — > C(e,T) at en- 
ergy e. If e = the Poincare map is the identity, so 
C(0, T) = Ct- For adiabatic invariance we need to prove 
that lim e ^o dC/de = 0, so that the difference between 
C{e,T) and Ct is of higher order than e [l4(- Since the 
contour C(e, T) can be locally represented by a function 
p x (x,e), we need to prove that lim e ^o dp x (x,e)/de = 0. 

In order to prove this, it is convenient to decompose 
the map Ct — > C(e, T) into three separate stages, start- 
ing out as an electron (from Ct to C+), followed by 
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FIG. 1: Classical trajectory in an Andreev billiard. Par- 
ticles in a two-dimensional electron gas are deflected by the 
potential V = [1 - (r/L) 2 ]V for r < L, V = for r > L. 
(The dotted circles are equipotentials.) There is specular re- 
flection at the boundaries with an insulator (thick solid lines) 
and Andreev reflection at the boundary with a superconduc- 
tor (dashed line). The trajectory follows the motion between 
two Andreev reflections of an electron near the Fermi energy 
Ef — 0.84 Vo- The Andreev reflected hole retraces this tra- 
jectory in opposite direction. 



FIG. 2: Poincare map for the Andreev billiard of Fig. Q 
Each dot represents a starting point of an electron trajectory, 
at position x (in units of L) along the interface y = and 
with tangential momentum p x (in units of ymVo). The inset 
shows the full surface of section, while the main plot is an 
enlargement of the central region. The drifting quasiperiodic 
motion follows contours of constant time T between Andreev 
reflections. The cross marks the starting point of the trajec- 
tory shown in the previous figure, having T = 18 (in units of 

y/mL a /Vo). 



Andreev reflection (C+ — > C_), and then concluded as 
a hole [from C_ to C(e,T)]. Andreev reflection intro- 
duces a discontinuity in p y but leaves p x unchanged, so 
C + = C-. The flow in phase space as electron (+) or 
hole (— ) at energy e is described by the action S±(q, e), 
such that p ± (q, e) = <9<S±/<9q gives the local dependence 
of (electron or hole) momentum p = (j> x ,p y ) on position 
q = (x,y). The derivative dS±/de — t±(q, e) is the time 
elapsed since the previous Andreev reflection. Since by 
construction t±(x, y = 0, £ = 0) = T is independent of 
the position x of the end of the trajectory, we find that 
lim e ^o dp x (x, y = 0, e)/de = 0, completing the proof. 

The drift (5x, 8p x ) of a point in the Poincare map is 
perpendicular to the vector (dT/dx, dT/dp x ). Using also 
that the map is area preserving, it follows that 

(Sx, S Px ) = ef(T)(dT/d Px , -dT/dx) + 0(e 2 ), (I) 

with a prefactor f(T) that is the same along the entire 
contour. 

The adiabatic invariance of isochronous contours may 
alternatively be obtained from the adiabatic invariance of 
the action integral I over the quasiperiodic motion from 
electron to hole and back to electron: 

/ = j pdq = e j j = 2eT. (2) 



Since e is a constant of the motion, adiabatic invariance 
of / implies adiabatic invariance of the time T between 
Andreev reflections. This is the way in which adiabatic 
invariance is usually proven in text books. Our proof 
explicitly takes into account the fact that phase space 
in the Andreev billiard consists of two sheets, joined in 
the constriction at the interface with the superconductor, 
with a discontinuity in the action on going from one sheet 
to the other. 

The contours of large T enclose a very small area. This 
will play a crucial role when we quantize the billiard, so 
let us estimate the area. It is convenient for this estimate 
to measure p x and x in units of the Fermi momentum pp 
and width W of the constriction to the superconductor. 
The highly elongated shape evident in Fig. [5] is a conse- 
quence of the exponential divergence in time of nearby 
trajectories, characteristic of chaotic dynamics. The rate 
of divergence is the Lyapunov exponent A. (We consider 
a fully chaotic phase space.) Since the Hamiltonian flow 
is area preserving, a stretching l + (t) = £ + (0)e xt of the 
dimension in one direction needs to be compensated by 
a squeezing l-{t) — £_(0)e~ At of the dimension in the 
other direction. The area A ~ £ + £_ is then time inde- 
pendent. Initially, £±(0) < 1. The constriction at the su- 
perconductor acts as a bottleneck, enforcing £±{T) < 1. 
These two inequalities imply £+(t) < e x ^~ T \ £_ < e~ At . 
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FIG. 3: Projection onto the x-y plane of the invariant torus 
with T = 18, representing the support of the electron compo- 
nent of the wave function. The flux tube has a large width 
near the superconductor, which is squeezed to an indistin- 
guishably small value after a few collisions with the bound- 
aries. 

The enclosed area, therefore, has upper bound 

A max ~ p F We~ XT ~ HNe~ XT , (3) 

where N ~ ppW/h 3> 1 is the number of channels in the 
point contact. 

We now continue with the quantization. The two in- 
variants e and T define a two-dimensional torus in the 
four-dimensional phase space. Quantization of this adi- 
abatically invariant torus proceeds following Einstein- 
Brillouin- Keller by quantizing the area 

<j) pdq = 2Trh(m + m = 0,1,2,... (4) 

enclosed by each of the two topologically independent 
contours on the torus. Eq. Q ensures that the wave- 
functions are single valued. (See Ref. |15| for a derivation 
in a two-sheeted phase space.) The integer v counts the 
number of caustics (Maslov index) and in our case should 
also include the number of Andreev reflections. 

The first contour follows the quasiperiodic orbit of Eq. 
©, leading to 

eT = (m+\)Trh, m = 0,1,2,... (5) 

The quantization condition is sufficient to determine 
the smoothed density of states p(e), using the classical 
probability distribution P(T) oc exp(-TNS/h) for 
the time between Andreev reflections. (We denote by 8 
the level spacing in the isolated billiard.) The density of 



states 

p{e) = N dT P{T) S(e - (m + ^U/T) (6) 

has no gap, but vanishes smoothly oc exp(— N8/4e) at 
energies below the Thouless energy N8. This "Bohr- 
Sommerfeld approximation" [12J has been quite success- 

fui mums, but it gives no information on the location 

of individual energy levels — nor can it be used to deter- 
mine the wave functions. 

To find these we need a second quantization condition, 
which is provided by the area <f T p x dx enclosed by the 
contours of constant T(x,p x ), 

j> Px dx = 2nh(n + v/4), n = 0,1,2,... (7) 

Eq. Q amounts to a quantization of the period T, which 
together with Eq. (0 leads to a quantization of e. For 
each T n there is a ladder of Andreev levels e nm = (m + 
\)-K%/T n . 

While the classical T can become arbitrarily large, the 
quantized T n has a cutoff. The cutoff follows from the 
maximal area (j3J enclosed by an isochronous contour. 
Since Eq. (JJJ requires j4 max > 2nh, we find that the 
longest quantized period is T = A _1 [lniV + 0(1)]. The 
lowest Andreev level associated with an adiabatically in- 
variant torus is therefore 

nh i:Ti\ 
£oo = W = 2frN- (8) 

The time scale To oc | ln?i| represents the Ehrenfest time 
of the Andreev billiard, which sets the scale for the exci- 
tation gap in the semiclassical limit [2(1 l2ll l22j | . 

We now turn from the energy levels to the wave func- 
tions. The wave function has electron and hole compo- 
nents i/)±(x, y), corresponding to the two sheets of phase 
space. By projecting the invariant torus in a single sheet 
onto the x-y plane we obtain the support of the elec- 
tron or hole wave function. This is shown in Fig. |3 for 
the same billiard presented in the previous figures. The 
curves are streamlines that follow the motion of individ- 
ual electrons, all sharing the same time T between An- 
dreev reflections. (A single one of these trajectories was 
shown in Fig. [I]) 

Together the streamlines form a flux tube that repre- 
sents the support of ip+ . The width SW of the flux tube is 
of order W at the constriction, but becomes much smaller 
in the interior of the billiard. Since 8W/W <£++£- < 
e A(t-T) _|_ e -A* ( witn o < t < T), we conclude that the 

flux tube is squeezed down to a width 

SW min ~ We- XT ' 2 . (9) 

The flux tube for the level eoo has a minimal width 
£Wmin — W/y/N. Particle conservation implies that 
| ■0+1 2 oc 1/8W, so that the squeezing of the flux tube 



4 



is associated with an increase of the electron density by 
a factor of y/N as one moves away from the constriction. 

Let us examine the range of validity of adiabatic quan- 
tization. The drift 5x, 8p x upon one iteration of the 
Poincare map should be small compared to W,pf- We 
estimate 



Sx 
W 



8p x 

PF 



nxN 



p -\(T -T n ) 

( m + h)^f — ■ (io) 



For low- lying levels (to ~ 1) the dimensionless drift is 
< 1 for T n < T Q . Even for T n = T Q one has Sx/W ~ 
1/lnJV < 1. 

Semiclassical methods allow to quantize only the tra- 
jectories with periods T < T . The part of phase space 
with longer periods can be quantized by random-matrix 
theory (RMT), according to which the excitation gap 
.Egap is the inverse of the mean time between Andreev 
reflections in that part of phase space 0, 0] : 



E — ^l 2 f) 

& gap — 7 It 



J™P(T)dT 

/" TP(T) dT T + 2irh/N6 



(11) 



Here 7 = — 1) is the golden ratio. This for- 

mula describes the crossover from E gajp = j 5 / 2 h/To = 
j 5 / 2 h\/\nN to E gap = 7 5 / 2 7V(5/27r at Am A ~ hX/5. It 
requires hX/NS ^> 1 (mean dwell time large compared to 
the Lyapunov time). The semiclassical (large- A) limit 
of Eq. (jTTJl . limAr^oo i? gap = 0.30 h/T is a factor of 5 



below the lowest adiabatic level, £00 = I.QU/Tq, so that 
indeed the energy range near the gap is not accessible by 
adiabatic quantization psT ]. 

Up to now we considered 2-dimensional Andreev bil- 
liards. Adiabatic quantization may equally well be ap- 
plied to 3-dimensional systems, with the area enclosed 
by an isochronous contour as the second adiabatic in- 
variant. For a fully chaotic phase space with two Lya- 
punov exponents Ai , A2 , the longest quantized period is 
To = h(Xi + A 2 ) _1 In A. We expect interesting quantum 
size effects on the classical localization of Andreev levels 
discovered in Ref. 0, which should be measurable in a 
thin metal film on a superconducting substrate. 

One important challenge for future research is to test 
the adiabatic quantization of Andreev levels numerically, 
by solving the Bogoliubov-De Gennes equation on a com- 
puter. The characteristic signature of the adiabatic in- 
variant that we have discovered, a narrow region of en- 
hanced intensity in a chaotic region that is squeezed 
as one moves away from the superconductor, should be 
readily observable and distinguishable from other fea- 
tures that are unrelated to the presence of the super- 
conductor, such as scars of unstable periodic orbits [lOj. 
Experimentally these regions might be observable using a 
scanning tunnelling probe, which provides an energy and 
spatially resolved measurement of the electron density. 
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